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几何 距 优化 质心 结合 隶属 度 约 束 RFCM 的 脑 MRI 图 像 分 割 算法 
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摘 要 : 针对 现 有 图 像 分 割 算 法 聚 类 复杂 以 及 分 割 精 度 不 够 高 的 问题 ， 提 出 了 基于 几何 距 优化 质心 和 粗糙 模糊 C- 均 值 
(RFCM) 相 结合 的 医学 图 像 聚 类 分 割 算法 。 首 先 建立 软 集 表示 的 像素 集 ， 并 计算 每 个 像素 与 质心 之 间 的 距离 ， 然 后 基于 
像素 和 质心 之 间 的 最 小 距离 ， 将 像素 分 组 到 聚 类 中 。 为 了 将 软 集 应 用 到 粗糙 模糊 C- 均 值 中 ， 定 义 了 一 个 模糊 软 集 ， 
一 步 将 输入 图 像 转换 为 二 值 图 像 ， 通 过 计算 连通 区 域 的 几何 距 选择 适当 的 质心 。 最 后 利用 这 些 新 的 质心 计算 更 新 像素 
的 隶属 度 值 ， 从 而 完成 模糊 聚 类 划分 。 在 Allen Brain Atlas 等 三 个 医学 数据 库 上 评估 了 所 提出 混合 算法 的 性 能 ， 获 得 的 
Jaccards 系数 和 分 割 精度 (SA) 都 优 于 几 种 对 比 算法 。 实 验证 明 ， 提 出 的 聚 类 分 割 算 法 具有 良好 的 性 
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RFCM magnetic resonance brain image segmentation algorithm based on geometric distance 
Co optimization of centroid and membership constraints 
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Abstract: Aiming at the problem that the existing image segmentation algorithm is complex and the segmentation precision is 


not high enough, proposed a medical image clustering segmentation algorithm based on geometric distance optimization 


pe centroid and rough fuzzy C-means (RFCM) . First set up the set of pixels represented by the soft set and calculate the distance 


between each pixel and the centroid. The pixels are then grouped into clusters based on the minimum distance between the 


pixel and the centroid. In order to apply the soft set to the coarse fuzzy C-means, a fuzzy soft set is defined. The input image is 
further converted into a binary image, and an appropriate centroid is selected by calculating the geometric distance of the 
connected region. Finally, using these new centroids to calculate the membership value of the updated pixel, the fuzzy 
clustering is completed. The performance of the proposed hybrid algorithm was evaluated on three medical databases, such as 
Allen Brain Atlas. The obtained Jaccards coefficient and segmentation accuracy (SA) were better than several comparison 
algorithms. Experiments show that the proposed clustering segmentation algorithm has good performance. 
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没有 考虑 空间 条 件 和 边界 条 件 ， 因 此 对 噪声 和 不 充分 数据 很 敏 

图 像 分 割 是 图 像 分 析 和 处 理 中 一 项 关键 的 功能 ， 在 描绘 解 ” 感 。 为 了 考虑 空间 条 件 和 边界 区 域 ， 一 些 研 究 开始 把 软 集合 理 
训 学 、 医 学 图 像 诊断 、 组 织 液 提取 、 病 理学 局 部 化 、 治 疗 计 划 、 ” 论 纳 入 到 KM 和 FCM 框架 中 进行 图 像 分 割 申 ,还 有 一 些 文献 中 
功能 图 像 数 据 的 部 分 容积 校正 以 及 计算 机 集成 外 科 手 术 等 方面 。 提出 了 混合 算法 : 粗 烟 K 均值 时 、 粗 糙 模 糊 可 能 性 C-35 (8 91. 
有 着 巨大 的 应 用 上。 模糊 集 、 粗 糙 集 是 处 理 医学 图 像 不 确定 性 ”阴影 C- 均 值 风 等。 文献 [12] 提 出 了 一 个 粗糙 集 的 算法 ， 使 用 
和 模糊 性 的 有 效 工 具 ， 被 广泛 应 用 于 医学 图 像 分 割 内 。 软 集合 。 个 立 值 而 不 是 负 域 ， 来 找到 下 近似 和 上 近似 。 然 而 ， 由 于 像素 
里 论 是 处 理 不 确定 性 和 模糊 性 的 一 种 新 的 数学 方法 口 绝对 属于 下 近似 和 上 近似 ， 或 者 两 者 都 有 ， 因 此 聚 类 时 不 会 消 
经 典 的 -均值 聚 类 算法 或 者 硬 C- 均 值 聚 类 算法 都 是 基于  ” 除 像 素 。 文 献 [13] 提 出 了 软 粗糙 集 ， 用 以 定义 不 包含 负 域 的 下 
清晰 集 G， 其 中 每 一 个 像素 只 属于 一 个 聚 类 。 它 的 模糊 版 本 : 近似 和 上 近似 。 但 是 这 种 技术 使 用 聚 类 原型 的 随机 初始 化 ， 当 
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选择 错误 时 会 产生 不 准确 的 结果 。 文 献 [14] 提 出 了 一 种 计算 模 
糊 软 集 的 下 近似 和 上 近似 的 方法 ， 以 有 效 地 定义 边界 。 然 而 ， 
这 些 定义 是 针对 清晰 值 ， 并 不 能 处 理 边界 像素 ， 粗 糙 区 域 的 确 
定 也 异常 复杂 。 这 几 种 基于 粗糙 集 的 算法 ， 都 涉及 复杂 的 指令 
来 计算 近似 值 ， 并 且 基 于 单一 六 值 将 像素 分 成 粗糙 区 域 。 试 图 
依据 聚 类 的 初始 质心 ， 使 用 像素 的 最 小 和 最 大 隶属 度 山 来 计算 
闵 值 。 因 此 ， 如 果 选 择 了 不 恰当 的 初始 质心 ， 会 导致 不 准确 的 
结果 。 此 外 ， 在 文献 [15] 中 提出 了 将 直觉 模糊 集 和 粗糙 集 相 结 
合 的 混合 算法 。 这 些 算法 使 用 了 三 元 隶属 度 ， 表 示 从 属于 一 个 
聚 类 、 不 属于 一 个 聚 类 、 不 确定 是 否 属 于 一 个 聚 类 。 使 用 三 元 
隶属 度 来 计算 上 下 近似 和 区 域 。 因 此 ， 这 些 算 法 涉及 的 计算 比 
粗糙 集 更 复杂 。 

本 文 提 出 了 使 用 软 集 表 示 图 像 的 观点 ， 并 用 软 模糊 粗糙 近 
似 定 义 了 图 像 粗 糙 区 域 的 下 近似 和 上 近似 。 与 粗糙 集 算 法 不 后 
这 种 聚 类 算法 避免 了 阀 值 .权重 参数 以 及 复杂 上 下 近似 的 使 用 。 


该 算法 不 使 用 负 区 域 ， 因 此 所 有 的 像素 都 可 以 参与 聚 类 。 软 集 
利用 相似 性 度量 来 量化 不 确定 性 ， 这 个 不 确定 性 跟前 一 个 聚 类 
和 当前 聚 类 的 相似 程度 有 关 。 
NE 粗糙 集 理论 

粗粮 集 对 不 精确 概念 的 描述 方法 是 通过 下 近似 和 上 近似 概 


念 来 描述 。 使 | 


粗粮 集 下 近似 和 上 近似 来 近似 不 确定 数据 。 设 


U 表示 全 集 ，R 是 等 价 关系 , U IREN DENA fa) 


的 集合 ， 它 们 在 X 中 形成 分 区 。(U,R) 是 近似 空间 。 对 子 集 
F 近 似 和 上 近似 可 以 表示 为 
Ra)= LU x 
XeX 


R(x,) = U X, 


x (X *ó 


(1) 


RO) 表示 对 象 X 属于 和 时 ，X 的 下 近似 空间 。R(x) 表 


示 对 象 可 能 属于 义 时，X 的 上 近似 空 
类 为 三 个 不 同 的 区 域 : 


Ep]. X 的 近似 空间 被 分 


正 域 (x,)= R(x) 
边界 域 (x,)= R(X) - RGs) Q) 
PR) =U 一 R(%) 


一 些 研 究 中 利用 粗粮 集 理论 ， 提 出 了 一 些 混合 算法 ， 如 : 
粗糙 K- 均 值 算 法 和 粗糙 模糊 C- 均 值 算 法 。 如 果 一 个 像素 属于 
所 有 的 聚 类 ， 则 其 隶属 度 较 少 。 把 该 像素 包括 到 聚 类 中 ， 会 产 
生 噪 声 ， 从 各 聚 类 中 去 除 该 像素 ， 会 产生 不 足 的 数据 。 聚 类 中 
的 这 些 不 完整 性 、 不 确定 性 和 模糊 性 可 以 使 用 粗粮 集 的 下 近似 
和 上 近似 来 定义 。 聚 类 的 下 近似 是 绝对 属于 聚 类 的 像素 集合 ， 
而 上 近似 是 可 能 属于 聚 类 的 像素 集合 。 给 予 边界 区 域 (上 下 近 
似 之 间 ) 中 的 像素 第 二 次 机 会 ， 以 便 在 下 一 次 迭代 中 它们 可 以 
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移动 到 合适 的 聚 类 。 因 此 ， 粗 糙 集 可 以 减少 聚 类 错误 。 
2 ”粗糙 模糊 C- 均 值 聚 类 


模糊 C- 均 值 聚 类 算法 是 从 硬 C- 均 值 聚 类 算法 发 展 而 来 , 它 
采用 和 迭代 法 优化 目标 函数 来 获取 对 数据 集 的 模糊 分 类 ， 算 法 具 
有 很 好 的 收敛 性 ， 是 实际 应 用 最 为 广泛 的 一 种 模糊 聚 类 算法 。 

文献 [16] 提 出 一 种 用 于 磁 共 振 脑 图 像 的 模糊 C- 均 值 聚 类 方 
法 ， 其 偏 置 场 模型 的 定义 如 方程 (3) 所 示 。 以 下 等 式 给 出 了 分 
割 磁 共 振 脑 图 像 的 目标 函数 ; 


J Gpem pa A ep 
Xem > (url 


xeB(C;) 


(3) 
2 
e 


其 中 :2 隶属 度 值 是 用 下 面 的 带 偏 置 场 的 等 式 计算 0 。 


| 
Ix; - 5C; (4) 
a i^ He 


MUS PAN BUE f, RUE, HER PCC) mih st B(C;) . 
Xin ^ BUE UTF: 


1 n 
f, m " 2. d i, 
"AS 
2 n^ max 


其 中 : d 是 像素 与 质心 的 最 小 距离 ，Qd 
大 距离 。 


(5) 


是 像素 与 质心 的 最 


P(Cjs Fid, < 
x 74B(Cjs Xin <d; <t, (6) 
N(C,; 其 他 
数 被 估计 为 正 域 与 整个 区 域 的 百分比 ， 


其 中 :每 个 聚 类 的 权重 参 
下 式 给 出 ; 


rc) 
=] 


W, = (7) 


=argmax (u,:i 


距离 dj 由 下 式 给 出 : 


d; = arg min(]x, -cD (8) 


E 交 多 项 式 的 线性 组 合 来 逐个 像素 地 近似 偏 置 场 b 。 


[ea 
E 


M 
b, = w,g,G) 2W'G() (9) 


其 中 :M =(D+D(D+D 是 度 为 DD 的 多 项 式 g, G) 的 数量 ， 而 
WOO 是 基 函 数 集合 G(i) 的 线性 组 合 。 使 用 能 量 函 数 估 计 和 
校正 区 域 Q 中 像素 i 的 偏 置 场 ， 如 下 式 所 示 : 


- AH 
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的 值 为 1.5> 1 <3， 能 够 获得 很 好 的 结果 ， 通 常 取 2。W 和 
W, =1 一 W 是 权重 因子 ， 影 响 下 近似 和 边界 区 域 的 相对 重要 


h 


F= SIl -w'Gac,[w a0) 
j ieQ 


能 量 相 对 于 人 最 小 化 ;可 以 通过 求解 DF / DW =0% 性 。 每 个 聚 类 都 与 下 近似 R(C ). 上 近 ri R(C, ) 和 边界 区 域 
完成 这 一 步 ， 得 到 : 


Do +24 ap — KC)-RC)-R(C,) 相关 联 。X 中 的 所 有 元 素 都 是 通过 使 
其 中 : A 是 由 下 式 给 出 的 多 维 列 向 量 ; 用 一 个 阔 值 来 分 组 到 下 近似 区 域 及 (C)) 或 边界 区 域 召 (C)) 或 
A=} GOGA «c (12) » 
气 上 近似 区 域 R(C,) P, LRL HH oo ox, 到 聚 类 质心 C) 的 
Ves xGü) «c, (13) — " l 
气 最 高 隶属 度 us 来 计算 的 。 粗 粮 模 糊 C- 均 值 使 用 式 〈17) 中 定 
HFE = ZR AJ — ^h £k TE 7. F: = 。 通 ? / 
4E DFIDW 一 0 表 不 为 一 个 线性 方程 AW-V 38 。 广 的 闭 值 ， 将 像素 置 于 下 近似 或 上 近似 中 。 


过 解 W = ACV WET MEHA b, =W GG) 。 在 这 种 方 

法 中 ， 由 闵 值 t 和 来 定义 粗糙 区 域 ， 距 离 大 于 的 像素 落 ris nti (17) 
arg max(u..) 

入 负 域 ， 因 此 从 聚 类 中 消除 该 像素 。 


MRTA, WAxeR(C)JERx eR), BW 


D 


3 ”提出 的 方法 


d 


g 
| 


Bes x ERC) EX e R(C,) 。 如 果 像 素 满足 阔 值 条 件 ， 那 么 
— 31.1 粗糙 模糊 C- 均 值 聚 类 "p : ET ES n 
e E PHEHENDEGE ATE BOEGE. 为 了 把 数据 。 它们 将 被 放置 在 下 近似 中 ， 否 则 它们 将 被 放置 在 边界 区 域 中 或 
e Md MMC MN 者 被 给 予 第 二 次 机 会 。 然 而 ， 由 于 缺乏 参数 化 工具 ， 使 用 粗糙 
O Pot AMEME CCC], WR FCM 集 计算 近似 值 是 比较 复杂 的 。 而且， 六 信和 权重 参数 是 用 户 定 
km 算法 的 抗 噪 性 ， 我 们 在 FCM 的 目标 函数 中 对 隶属 度 引 入 空间 — nfi. 
— 约束 关系 ， 新 的 目标 函数 为 3.1.2 模糊 软 集 
2 j- (U, v) Sur (i, j)d? (i. j)+ 软 集 是 将 参数 映射 到 全 集 。 令 U 为 全 集 ， A 是 参数 集 ， 

k=l ij 

a) EBACE. 那么 U 上 的 软 集 由 (了 ,A) 给 出 ,其 中 下 是 一 
ZIP) MM 
"AE 个 映射 ，F: AS P(U). PUERU l4, AU 的 


EP: N, 为 像素 点 的 邻 域 区域 ， 上 式 的 第 二 项 对 (站 zog ”所 有 了 集 “。 把 软 集 应 用 于 图 像 分 别 首 先 把 观察 图 像 中 的 像素 


ai Ui (i j) ts3tsuse mu, (I, j) 的 差异 进行 惩罚 。 为 。 集 表示 为 U ， U-X ={z|5 是 图 像 中 第 ;个 像素 的 值 } 。 候 


求解 隶属 度 函数 ， 建 立 Lagrange 函数 ， 可 得 定 观察 图 像 X 是 由 真实 信号 通过 偏 置 场 b 和 零 均值 高 斯 噪声 
C 
F = (i 7) dr (5 3) * 4. [站 得 到 的 。 然 后 设 (五, A) 是 将 要 构建 的 软 集 ,包含 属于 聚 类 的 像 
d (15) 
素 。 最 后 把 参数 集 作为 像素 所 属 聚 类 C{i —1,2,3, .. kic 
PP) XA 


假设 图 像 中 的 像素 集 由 UU = Dux, XA Hs. SURE 


Ui 表示 第 i 个 像素 对 第 j 个 模糊 聚 类 的 隶属 程度 ， 由 式 
A 为 像素 所 属 聚 类 {C,C,,C;}。 计算 每 个 像素 与 质心 之 间 的 距 


(16) 给 出 。 bk 
1 离 。 基 于 像素 和 质心 之 间 的 最 小 距离 ， 按 如 下 方法 将 像素 分 组 
u; = — 到 聚 类 中 。 
EI x —C, |" (16) 
iEcel gl FG) - bs] 
e s F(C,)= {X1 X3, X4, Xs} (18) 


F(C;) —[x4,X4,X;, X4] 


C ERX j 的 平均 值 ，1 是 用 户 定 义 的 模糊 值 。 如 果 模 


上 面 例子 的 软 集 表示 为 


糊 值 m 接近 于 1, 隶属 度 函数 变 得 更 加 清晰 明了 , 并 且 如 果 m 
增加 , 隶属 度 函 数 变 得 越 来 越 模糊 不 清 。 对 于 大 部 分 数据 ，m 
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Ci = {%, Xe} 
(F, A) 4C, ={X,X,, X4: X5} (19) 
C, = {%, X4, Xs, X4] 
表 1 为 基于 软 集 表示 的 像素 集 。 WR, e F(C;), WA 
其 值 为 1， 否则 就 是 0。 
dl 基于 软 集 表示 的 像素 集 
U 1 C, C, 
X 0 1 1 
X 1 0 0 
X 0 1 1 
Xa 0 1 0 
Xs 1 0 1 
Xe 1 1 1 
T 为 了 将 软 集 应 用 到 模糊 值 中 ， 定 义 了 一 个 模糊 软 集 ; 模糊 
co KE (F.A) 是 U 的 模糊 值 子 集 的 参数 秘 。 DRUG, (F.A) 可 以 
~ 
E REH (FP, A) -LFCO | Ge AY. {FORRU 的 模糊 子 集 ， 
z EJ E BORA es FU) = {( Ar )} 。 因 此 ，X 属 于 
Co U ， 其 隶属 度 Wj(X) 为 参数 。 因 此 ， 通 过 考虑 每 个 像素 对 聚 类 
L 质心 的 模糊 隶属 度 ， 以 软 模糊 集 表示 图 像 。 模 糊 软 集 是 参数 对 
CO 。 具有 模糊 值 的 全 集 映射 。 因 此 ， 为 了 表示 模糊 软 集 ， 使 用 等 
p (200 计算 像素 的 隶属 度 值 。 假设 像素 对 聚 类 的 隶属 度 值 如 下 : 
s F(C,) = ((x,,0.6),(x,,0.1) (x,,0.4),(x,.0.7)) 
人 F(C,) = ((3,.0.4),(2,,0.2). (35.0.5). (x,,0.1) Q0) 
= F(C,) = ((x,,0.4) (25.0.3) (3.0.5); (x,,0.2)] 
É m 
@O 5b, BUB: u= (C,.0.6) (Cs,0.6) (C3,0.6， 这 样 一 个 
0.6 隶属 度 的 像素 能 被 放置 在 聚 类 中 。 选择 0.6 的 值 是 为 了 包括 
在 聚 类 中 的 像素 至 少 有 0.6 的 隶属 度 值 。 如 果 它 小 于 0.5， 那 么 
聚 类 会 变 得 更 加 模糊 并 增加 噪声 。 关于 软 集 的 人 模糊 粗糙 的 下 
近似 和 上 近似 由 RA(X) 和 Ru) 定义 。 


RUDA a (ai G CO G)) (A (1 F0) G))) 


(21) 


RUDA a (FO) (vs ((F CO) a) 2 


和 (22) 


可 以 修改 成 如 下 


EU 


于 图 


把 U 作为 X，K 作为 模糊 集 ，U; EIFU), 方程 QD 
像 分 割 。 


RUA ec, (1-1) v (es (11) v G))) 


Q3) 


RUA, «c, (Gu) vae (us) uG))) Q4) 
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南 丽 丽 ， 等 : 几何 距 优 化 质心 结合 来 属 度 约束 RFCM 的 脑 MRI 图 像 分 割 算法 
使 用 等 式 (23) 和 24) 的 模糊 集 KH 计算 软 模糊 集 的 下 近 


似 和 上 近似 ， 下 近似 和 上 近似 的 表格 表示 在 表 2 中 给 出 。 


表 2 下 近似 值 和 上 近似 值 
U C, C, C, Ru(x) Ru(x) 
X, 0.5 0.6 04 04 0.6 
X; 0.1 0.2 0.2 0.7 0.1 
x, 04 0.6 0.5 0.4 0.4 
x; 0.6 0.1 0.2 0.5 0.2 
H 0.6 0.7 0.6 — = 
在 上 面 的 定义 中 , 没有 使 用 阔 值 或 权重 参数 来 计算 近似 值 ， 
而 是 使 用 了 简单 的 与 、 或 操作 。 使 用 这 些 近似 值 ， 就 可 以 通过 
以 下 等 式 确定 正 域 P(C, ) 和 边界 域 B(C ,)。 
(x) 
(25) 


P(C); 若 x, > Ryu 
Ru(x,) 


(BC): ZRU) x, < 
3.1.3 几何 距离 优化 质心 
如 果 一 个 像素 的 度数 在 Ru(x) 内 ， 那 么 把 该 像素 置 于 边界 
域 ， 否 则 ， 如 果 该 像素 具有 比 Ru) 更 大 的 值 ， 则 把 该 像素 置 
于 正 域 中 。 一 旦 获得 了 边界 域 和 正 域 ， 那 么 就 可 以 基于 正 域 和 
负 域 ， 使 用 下 面 公式 中 的 偏 置 场 b; 来 计算 新 的 质心 。 


b3 xb, > (u, y xb; 
xeMc) — ) nS 
C,- —€—À———— (6) 
Ei b; bi (wu, ) b? 
xeP(c;) xeb(c;) 
利用 这 些 新 的 质心 ， 计 算 更 新 像素 的 隶属 度 值 。 由 于 聚 类 
算法 是 迭代 的 ， 基 于 软 集 的 聚 类 算法 使 用 前 一 次 迭代 和 当前 过 
代 形 成 的 聚 类 之 间 的 相似 度 进行 迭代 。 
初始 聚 类 质心 在 分 割 中 扮演 着 重要 的 角色 ， 在 初始 隶属 度 
函数 中 使 用 如 K- 均 值 聚 类 算法 或 者 硬 C 均值 聚 类 算法 这 样 的 
聚 类 算法 。 错 误 选 择 的 质心 可 能 导致 局 部 最 小 从 而 导致 不 准确 
的 结果 。 将 输入 图 像 转换 为 二 值 图 像 ， 并 选择 相应 的 灰 度 作为 


聚 类 质心 。 几 何 距 优化 质心 的 算法 步骤 如 下 : 
引 将 输入 图 像 转换 为 二 值 图 像 ; 
b) 通 过 连通 组 件 标记 算法 找到 所 


初始 村 


了 的 连通 区 域 ， 


分 别 标 


ii 


c) 对 每 个 连通 区 域 运用 计算 几何 距 算法 得 到 质心 ; 


M, =f fix yf( (x. y )dxdy 
不 同 颜色 绘制 出 连通 区 域 与 质心 ， 


Q7) 


进一步 输出 处 理 


rus 


d) 再 
后 图 像 ; 
软 模糊 粗糙 外 
到 的 。 软 粗糙 集 使 ) 
近似 和 上 近似 。 每 
来 定义 的 。 


前 过 杉 入 软 集 、 粗 糙 集 和 模糊 集 的 概念 得 
表示 归属 程度 的 软 模糊 集 来 定义 聚 类 的 下 
近似 值 都 是 用 模糊 值 而 不 是 清晰 值 
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在 所 有 大 脑 图 像 的 实验 设置 中 ， 取 聚 类 数量 上 =4〈 关 于 

GM, WM, CSF 和 背景 ) 。 使 用 与 实际 分 割 边 界 的 定量 比较 来 

4.1 方法 实施 计算 分 割 结果 ， 使 用 Jaccards 系数 (Jaccard Coefficient, JC) 和 分 
在 三 个 医学 数据 库 上 评估 了 所 提出 混合 算法 的 性 能 : G) 割 精 度 (Segmentation Accuracy, SA) 对 分 割 结果 进行 评估 。 两 个 
脑 图 像 数据 库 Brain Perfusion Database!?; Gi) 艾 伦 大 脑 图 谱 。 ”数据 集 之 间 的 JC 是 两 个 数据 集 共 有 的 属性 数量 除 以 所 有 属性 
Allen Brain Atlas"?!; (iii) BRATS 数据 库 Po0。 这 些 数据 库 中 选 ”数量 的 结果 。 


4 ”实验 结果 


取 的 供 仿真 所 用 图 像 的 数量 以 及 大 小 如 表 3 所 示 。 这 三 个 数据 „u ANB 

i Jacca1rd 系 数 = 一 一 (29) 
库 都 是 经 典 的 脑 图 像 数据 库 ， 图 像 的 数量 以 及 大 小 都 适合 本 文 AUB 
方法 进行 仿真 实验 。 把 这 些 图 像 转 换 成 MATLAB 可 读 的 格式 ， 其 中 : A 是 本 文 方法 获得 的 分 割 图 像 ，B 是 真实 图 像 。 
作为 预 处 理 步 又， 在 进行 分 割 之 前 ， 使 用 MIPAV〔 医 学 图 像 处 Jaccard 系数 在 70% 以 上 ， 这 意味 着 分 割 结果 是 好 的 。SA 是 用 


里 、 分 析 和 可 视 化 ) 软件 将 非 脑 部 图 像 从 图 像 中 移 除 。 所 提出 来 比较 分 割 图 像 和 真实 图 像 之 间 相 似 性 的 另 一 个 常用 指标 。 为 
的 算法 在 4.2GHz Intel I7 7700K 处 理 器 上 的 MATLAB 2012a 中 了 计算 分 割 精度 ， 对 每 个 分 割 结 果 计算 四 个 参数 。 
实现 。 图 1 为 艾 伦 大 脑 图 谱 中 的 部 分 图 像 。 a) 真 正 (TP): 真实 图 像 中 正确 像素 被 正确 检测 为 分 割 像素 。 
b) 真 负 (TN): 真实 图 像 中 错误 像素 被 正确 检测 为 分 割 像素 。 
cMBAE(FP): 分 割 区 域内 未 发 现 真 实 图 像 中 的 正确 像素 。 
EHEN): 分 割 区 域 中 不 存在 真实 图 像 中 的 错误 像素 。 
| TP« FN 
OTP& FP+TN + FN 
第 一 个 实验 在 脑 图 像 数据 库 Brain Perfusion Database 数据 
库 上 进行 ， 该 数据 库 包 含 TI 加 权 仿 真 大 脑 10 5 E — 465 
向 视图 ， 片 数 为 90， 厚 度 为 Imm， 不 同 强度 不 均匀 性 ANU = 
0 和 INU=20) 以 及 噪声 级 别 (0%、1%、3%、5% 和 9%)。 
于 磁 共 振 脑 图 像 数据 的 典型 噪声 级 别 在 3% 左 右 ， 因 此 选择 


(30) 


图 1 艾 伦 大 脑 图 谱 中 的 部 分 图 像 含有 3% 噪 声 和 0%INU 的 合成 脑 磁 共振 图 像 。 为 了 验证 实验 ， 
表 3 仿真 所 用 数据 集 还 选择 了 其 他 噪声 级 别 为 9% 的 噪声 图 像 。 
维度 数据 集 数 图 像 大 小 4.2 性 能 评估 
2D Brain Perfusion Database 10 194*237 许多 聚 类 算法 的 性 能 受 初始 聚 类 质心 的 影响 。 本 文 提出 的 
2D Allen Brain Atlas 10 256*128 直方 图 优化 质心 的 方法 解决 了 关于 随机 质心 的 问题 。 基 于 直方 
2D BRATS 10 变量 *216 图 的 优化 质心 在 寻找 接近 真实 质心 的 起 始点 方面 实现 了 高 效 
磁 共 振 脑 图 像 经 常 分 为 白质 (WM) 、 灰 质 (GM) po — 率 ， 从 而 得 到 可 靠 的 脑 图 像 分 割 结果 。 表 4 给 出 了 算法 的 SA 


液 《CSF) 。 有 效 提取 这 些 区 域 对 于 定量 分 析 是 非常 重要 的 ， 结果 ， 表 5 中 给 出 了 算法 的 JC 值 。 结 果 表 明 ， 即 使 存在 噪声 ， 
因此 ， 产 生 了 分 制 MR 脑 图 像 的 几 种 方法 。 磁 共振 脑 图 像 如 图 ”提出 的 新 算法 也 表现 得 非常 好 。 

2 所 示 。 文献 [12] 方 法 虽然 花费 很 少 的 时 间 进 行 评 估 ， 但 对 存在 于 
边界 区 域 〈 不 可 辨认 ) 的 像素 却 不 起 作用 ， 因 此 增加 了 聚 类 错 
WR, 降低 了 SA 和 JC 值 。 虽然 通过 定义 隶属 度 ，K 均值 聚 类 算 
或 者 硬 C- 均 值 聚 类 算法 方法 能 够 来 处 理 模糊 性 ,但 对 不 可 分 
的 情况 却 不 能 恰当 地 聚 类 。 因 此 , 一 个 组 织 的 SA 和 JC 值 
0， 男 一 个 组 织 的 SA 和 JC 值 就 会 减 小 。 通过 提供 充当 缓冲 
域 的 边界 区 域 ,粗糙 集 可 以 有 效 地 处 理 聚 类 中 的 不 可 分 辨 像素 。 
在 分 配 到 聚 类 之 前 ， 给 予 缓冲 区 中 的 像素 第 二 次 机 会 ， 以 减少 
聚 类 错误 ， 获 得 良好 的 SA 和 JC 结果 ,如 表 4 表 5 以 及 图 3 所 
示 。 而 文献 [13-15] 方 法 基于 单一 阔 值 将 像素 分 成 粗糙 区 域 。 试 
图 依据 聚 类 的 初始 质心 ， 使 用 像素 的 最 小 和 最 大 隶属 度 来 计算 
阔 值 。 因 此 ， 如 果 选 择 了 不 恰当 的 初始 质心 ， 会 导致 不 准确 的 
图 2 磁 共 振 脑 图 像 结果 。 


这 


t 


Ico ub o n 
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dA 提出 的 方法 与 不 同 分 割 算法 在 不 同 噪声 级 别 脑 图 像 上 仿真 的 性 能 比较 (分 割 精度 ) 


高 斯 噪声 组 织 类 型 文献 [12] 文献 [13] 文献 [14] 文献 [15] 提出 的 方法 
GM 0.8123 0.8256 0.8055 0.8427 0.8976 
3% CSF 0.8352 0.8692 0.8460 0.8583 0.8892 
WM 0.8276 0.8027 0.8289 0.8167 0.9289 
GM 0.7967 0.7866 0.6268 0.6822 0.7992 
6% CSF 0.6976 0.8055 0.7199 0.7620 0.8027 
WM 0.6678 0.6311 0.7278 0.7856 0.7945 
GM 0.8027 0.7289 0.6678 0.6822 0.7978 
9% CSF 0.6918 0.7389 0.7278 0.7620 0.7892 
WM 0.6899 0.7027 0.7901 0.7856 0.8267 
表 5 提出 的 方法 与 不 同 分 割 算法 在 不 同 噪声 级 别 脑 磁 共 振 图 像 上 仿真 的 性 能 比较 (Jaccards 系数 ) 
高 斯 噪声 组 织 类 型 文献 [12] 文献 [13] 文献 [14] 文献 [15] 提出 的 方法 
GM 0.7826 0.8018 0.7822 0.7266 0.8729 
3% CSF 0.8036 0.8026 0.7811 0.8177 0.8677 
WM 0.8207 0.8237 0.8026 0.8027 0.8366 
GM 0.7826 0.7941 0.6622 0.6972 0.8027 
6% CSF 0.7027 0.7832 0.7598 0.7620 0.8062 
WM 0.6388 0.6432 0.7801 0.7893 0.7945 
GM 0.7372 0.7941 0.6748 0.6822 0.8072 
9% CSF 0.7026 0.7832 0.7528 0.7620 0.8053 
WM 0.6892 0.7027 0.7729 0.7836 0.8366 
! aoxENA ü 
09 5% URE 25 
0.8 
0.7 a 全 20 
E | | | | | | | 
É dM 
0.3 
0.2 
P 文献 D2] XANI 。 文献 04] XAS “提出 的 方法 
文献 [12] 。 文献 [13] — 文献 [14] — 文献 [15] ”提出 的 方法 到 4 不 同 分 审 算 法 需要 的 执行 时 间 比 较 
图 3 不 同 分 割 算法 在 不 同 噪声 级 别 脑 图 像 上 仿真 
的 性 能 比较 (分 割 精 度 ) 5 in 
软 集 使 用 参数 化 工具 有 助 于 有 效 设计 上 近似 和 下 近似 。 软 本 文 提 出 了 基于 粗糙 模糊 C- 均 值 隶属 度 约 束 的 磁 共 振 脑 
集 也 使 用 相似 性 系数 ， 将 一 个 数量 分 配给 代表 相等 程度 的 一 对 “图像 聚 类 分 割 算 法 。 将 软 集 应 用 于 模糊 集 和 粗糙 集 ， 提 出 了 混 
软 集 。 与 许多 聚 类 算法 类 似 ， 基 于 软 集 的 算法 也 是 迭代 的 ， 合算 法 及 其 在 医学 图 像 分 割 中 的 应 用 。 提 出 的 算法 具有 许多 优 
且 由 于 软 集 使 用 参数 化 工具 和 相似 性 系数 ， 因 此 可 以 非常 快速 ” 点 : 利用 参数 化 工具 定义 上 近似 和 下 近似 来 处 理 不 确定 性 ， 减 
地 达到 收 化 。 为 了 处 理 不 可 分 辨 的 像素 并 减少 时 间 ， 粗 糙 集 与 。” 少 聚 类 错误 ， 使 用 软 集 相 似 性 系数 达到 快速 收敛 。 该 算法 产 昌 


软 集 一 起 使 用 被 证 明 是 更 好 的 解决 方案 。 对 各 个 算法 的 执行 时 精确 的 结果 。 根 据 综合 数据 和 仿真 评估 所 提出 算法 的 性 能 ， 所 
间 进 行 比较 ， 结 果 如 图 4 所 示 。 提出 的 方法 提供 了 准确 的 大 脑 图 像 分 割 。 未 来 希望 能 够 对 提出 
于 使 用 参数 化 和 相似 性 系数 ， 所 提出 的 方法 避免 了 闵 值 的 方法 进行 进一步 完善 。 
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